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ABSTRACT 

Evidence suggesting an observable magnetic interaction between a star and 
its hot Jupiter {Porb < 7 days, a < 0.1 AU, Mpsini > 0.2 Mj) appears as a cyclic 
variation of stellar activity synchronized to the planet's orbit. In this study, 
we monitored the chromospheric activity of 7 stars with hot Jupiters using new 
high-resolution echelle spectra collected with ESPaDOnS over a few nights in 



^ Based on observations collected at the Canada- France-Hawaii Telescope operated by the National Re- 
search Council of Canada, the Centre National de la Recherche Scientifique of France, and the University of 
Hawaii. 
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2005 and 2006 from the CFHT. We searched for variabihty in several stellar 
activity indicators (Ca II H A3968, K A3933, the Ca II infrared triplet (IRT) 
A8662 line, Ha A6563 and He I A5876). HD 179949 has been observed almost 
every year since 2001. Synchronicity of the Ca II H & K emission with the orbit 
is clearly seen in four out of six epochs, while rotational modulation with Prot=7 
days is apparent in the other two seasons. We observe a similar phenomenon 
on V And, which displays rotational modulation (Proi= 12 days) in September 
2005, in 2002 and 2003 variations appear to correlate with the planet's orbital 
period. This on/off nature of star-planet interaction (SPI) in the two systems is 
likely a function of the changing stellar magnetic field structure throughout its 
activity cycle. Variability in the transiting system HD 189733 is likely associated 
with an active region rotating with the star, however, the flaring in excess of 
the rotational modulation may be associated with its hot Jupiter. As for HD 
179949, the peak variability as measured by the mean absolute deviation (MAD) 
for both HD 189733 and r Boo leads the sub-planetary longitude by ~ 70°. The 
tentative correlation between this activity and the ratio of MpSini to the planet's 
rotation period, a quantity proportional to the hot Jupiter's magnetic moment, 
first presented in Shkolnik et al. 2005 remains viable. This work furthers the 
characterization of SPI, improving its potential as a probe of extrasolar planetary 
magnetic fields. 

Subject headings: stars: late-type, activity, chromospheres, planetary systems, 
radiation mechanism: non-thermal, stars: individual: r Boo, HD 179949, HD 
209458, HD 189733, HD 217107, HD 149143 

1. Introduction 

Observations and theory demonstrate that star-planet interaction (SPI) is a complex, 
yet potentially very informative probe of extrasolar planetary magnetic fields. In Shkolnik 
et al. (2003, 2005a), we reported on planet-induced chromospheric activity on two stars, 
HD 179949 and v And apparent from the night-to-night modulation of the Ca II H & K 
chromospheric emission phased with the hot Jupiter's orbit. The modulation was indicative 
of a magnetic rather than tidal interaction (Cuntz et al. 2000), such that the period of the 
observed stellar activity correlated with the planet's orbital period Pgrb, rather than Porb/'^- 
Ample observational evidence of tidal and magnetic interactions exists in the exaggerated 
case of the RS Canum Venaticorum (RS CVn) stars, which are tightly-orbiting binary sys- 
tems consisting of two chromospherically active late-type stars (e.g. Glebocki et al. 1986, 
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Catalano et al. 1996, Shkolnik et al. 2005b ). 

Although efforts to observe variable radio emission from the stars with hot Jupiters have 
not yet been successful (e.g. Lazio & Farrell 2007, George & Stevens 2007), there have been 
several additional observations that support the existence of SPI. Photometric observations 
by the MOST space telescope of several hot- Jupiter systems, including HD 179949 and 

T Boo, suggest that stellar surface activity in the form of active spots may be induced by the 
giant planet (Walker et al. 2006, 2007). Also, Saar et al. (2006) recently reported a possible 
detection of planet-induced X-ray emission from the HD 179949 system corresponding to 
~30% increase in X-ray flux over quiescent levels coincident with the phase of the Ca 11 
enhancements at 0or6~O.8. A statistical analysis by Kashyap et al. (2006) suggests that 
the X-ray flux from stars with hot Jupiters is on average > 3 times greater than stars with 
planets at larger orbital distances, presenting further evidence that close-in giant planets 
have a measurable effect on the activity of the parent star. 

One scenario of magnetospheric interaction proposes that the planet induces reconnec- 
tion events as it travels through the large stellar magnetic loops (Cuntz et al. 2000, Ip et 
al. 2004), implying that the resulting activity should depend on the star's magnetic field, the 
planet's magnetic field and the orbital distance with respect to the Alfven radius of the host 
star (~10 stellar radii). Novel research of the magnetic field topology of hot Jupiter host 
stars is underway (Catala et al. 2007, Moutou et al. 2007) using Zccman Doppler Imaging 
(ZDI), which hopes to contribute to a more detailed understanding of SPI. 

A detection of a magnetic field of a hot Jupiter would 1) provide a constraint on the 

rapid hydrodynamic escape of its atmosphere (Vidal-Madjar et al. 2003, 2004) which could 
affect the planet's structure and evolution, 2) present implications for the planet's internal 
structure, and 3) shed light on the mass-radius relationship of the known transiting planets 
(Pont et al. 2005, Bakos et al. 2006). Although the internal magnetic fields of hot Jupiters 
are expected to be weaker than Jupiter's due to probable tidal locking and slower spin 
rates (Sanchez-Lavega 2004, Griessmeier et al. 2004), Olson & Christensen (2006) calculated 
that the magnetic field of a planet with even a tenth of Jupiter's rotation rate would still 
have a strong dipole moment, when reasonably assuming that the convection is not highly 
modified by the rotation rate. Also, the fact that both hot and very hot Jupiters, such as 
HD 209458 b and OGLE-TR-56 b, are detected at all means that they must have strong 
enough magnetic fields to balance the extreme stellar irradiation and CME plasma pressure 
to prevent destructive atmospheric erosion (Khodachenko et al. 2007). 

We seek to probe hot Jupiter magnetic fields in order to understand their formation and 
evolution. SPI potentially offers an indirect way to detect, and with future modeling and 
observations, measure planetary magnetic fields. 
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It is reasonable to assume that any magnetic interaction would be greatest in the out- 
ermost layers of the star, namely the chromosphere, transition region and the corona due to 
their proximity to the planet, low density, and non-radiative heat sources. With the commis- 
sioning of CFHT's high-resolution echelle spectrograph, ESPaDoNS, we are able to include 
several stellar activity indicators in our analysis to observe the interaction as a function of 
atmospheric height in order to model the energy transfer and dissipation mechanisms of this 
phenomenon. ESPaDOnS' wavelength coverage allows simultaneous monitoring of the Ca II 
infrared triplet (IRT, lower chromosphere). Ha, Ca II H, K (middle chromosphere), and He 
I D3 (upper chromosphere). 

Our program stars have planets with orbital periods between 2.2 and 7.1 days, eccentric- 
ities ~ and semi-major axes < 0.08 AU. These systems offer the best chance of observing 
upper atmospheric heating. Of the seven systems we observed with ESPaDoNS, r Boo, 
HD 179949, v And and HD 209458 have been observed previously in our CFHT/Gecko 
campaign. The first results from 2001 and 2002 observations, including the first evidence 
of planet-induced magnetic heating of HD 179949, were published in Shkolnik et al. (2003, 
2005a). We later extended the experiment at the Very Large Telescope (VLT) to include 
five southern targets. The three new systems monitored in this study are HD 217107, HD 
149143 and HD 189733. The system parameters for the ESPaDOnS program stars are listed 
in Table 1 along with our standard 61 Vir. 

In this paper, we present new echelle spectra and compare with those of previous years, 
bringing to light a broader understanding of stellar activity, its cycles, and SPI. The details 
of our observations and data reduction are outlined in Section [2l In Section [3] we discuss 
our analysis and results of the Ca II K measurements including long-term, short-term and 
rotational modulation. Comparisons with other activity indicators are made in Section HI 
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Table 1. Stellar and Orbital Parameters 



Star 


SpT 


vsini 


Prot 


^ orb 




a* 


(K)b 




(MADK)*! 


He I EW 






(km s-i) 


(days) 


(days) 


(Mj) 


(AU) 


(A) 


(A) 


(A) 


(ml) 


r Boo 


F7 IV 


14.8±0.3 


3.2*^ 


3.31 


4.4 


0.046 


0.336 


0.184 


0.0019 


27 


HD 179949 


F8 V 


6.3±0.9 




3.09 


0.98 


0.045 


0.369 


0.186 


0.0022j 


17 


HD 209458 


GO V 


4.2±0.5 


16S 


3.53 


0.69^ 


0.045 


0.195 


0.078 


0.0009 


^3 


V And 


F7 V 


9.0±0.4 




4.618 


0.71 


0.059 


0.254 


0.091 


0.0016 




HD 189733 


Kl V 


2.92±0.22 


11.7' 


2.22 


1.15^ 


0.031 


1.337 


1.231 


0.0044j 


35 


HD 217107 


G8 IV 


9.0±0.4 


39^ 


7.13 


1.35 


0.075 


0.160 


0.075 


0.0007 


<4 


HD 149143 


GO IV 


3.9±1 


?? 


4.09 


1.36 


0.052'' 


0.342 


0.144 


0.0014 


10 


61 Vir 


05 V 


2.2™ 


33"^ 








0.182 


0.083 


0.0008 


<4 



^Published orbital solutions: r Boo - Butler et al. 1997, HD 179949 - Tinncy ct al. 2001, HD 209458 - Charbonnoau 
et al. 1999, v And - Butler et al. 1997, HD 189733 - Bouchy et al. 2005, HD 217107 - Fischer et al. 1998, HD 149143 
- Dc Silva ct al. 2006 

''Total integrated intensity of the mean normalized Ca 11 K core. These values are relative to the normalization 
points near 3930 and 3937 A at 1/3 of the pseudo-continuum at 3950 A. 

^We subtracted the photospheric emission from (K) in order to measure the mean integrated chromospheric emission 
(K') using data from Wright et al. (2004). (See text for more details.) 

•^Average integrated 'intensity' of the mean absolute deviation (MAD) of the K residuals, per observing run 

''Henry et al. 2000 

^This work and Wolf & Harmanec 2004 

SMazch ct al. 2000 
''Transiting system 
'CroU et al. 2007 

JValuc were corrected to remove geometric (rotational and/or planetary) modulation of an active region on the star. 
For HD 189733, the non-corrected value is 0.0098 and for HD 179949, 0.0063. 

'^Wright et al. 2004 
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Table 2. Observations 



Star 


U 


B 


V 


Exposures'^ 


S/N^ 


S/N^ 










t X n X N 


at 3950 A 


at 8710 A 


r Boo 


5.02 


4.98 


4.50 


120 X 10 X 9^ 


630 


1640 


HD 179949 


6.83 


6.76 


6.25 


660 X 5 X 7" 


440 


1270 


HD 209458 


8.38 


8.18 


7.65 


1800 X 4 X 4 


350 


960 


V And 


4.69 


4.63 


4.09 


180 X 120 X 4° 


2800 


7290 


HD 189733 




8.60 


7.67 


1800 X 4 X 4 


280 


1260 


HD 217107 


7.33 


6.90 


6.18 


600 X 5 X 3 


350 


1340 


HD 149143 




8.53 


7.90 


1800 X 4 X 4 


320 


1000 


61 Vir 


5.71 


5.45 


4.74 


120 X 7 X 5 


450 


1480 



H = Exposure time in seconds, n = number of exposures per night, N = number 
of nights during the June 2006 ESPaDOnS observing run 

''Typical nightly S/N per 0.022-A pixel 

'^Spectra from three of the nine nights were observed in ESPaDOnS' 'spectropo- 
larimetry' mode. 

^Four nights in June 2006 and three nights in September 2005 

''All these data were acquired in the 'spectropolarimetry' mode in September 
2005. The extremely high S/N was a requirement of the partner program to search 
for linear polarization. (Collier Cameron et al., in prep.) 
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2. The spectra 

The observations were made with the 3.6-m Canada-France-Hawaii Telescope (CFHT) 
on 7 nights in September 2005 and 9 nights in 2006 June. We used ESPaDOnS (Echehe Spec- 
troPolarimetric Device for the Observation of Stars), which is fiber fed from the Cassegrain 
to coude focus where the fiber image is projected onto a Bowen-Wahaven slicer at the 
spectrograph entrance. With a 79 gr/mm grating and a 2048 x 4608-pixel CCD detector, 
ESPaDOnS' 'star-only' mode records the full spectrum over 40 grating orders covering 3700 
to 10400 A at a spectral resolution R of ~80,000. The four nights of observations of v And 
(18, 19, 23, 24 September 2005) and three nights of the r Boo data (16, 17, 18 June 2006) 
were taken in ESPaDOnS' 'spectropolarimetry' mode at R of 68,000. 

The data were reduced using Libre Esprit, a fully automated reduction package provided 
for the instrument and described in detail by Donati et al. (1997, 2007lj). Each stellar 
exposure is bias-subtracted and flat-fielded for pixel-to-pixel sensitivity variations. After 
optimal extraction, the 1-D spectra are wavelength calibrated with several Th/Ar arcs taken 
throughout the night. Finally the spectra are divided by a flat-field response and then 
the continuum is normalized. Heliocentric velocity corrections are applied as well as small 
velocity corrections (< 100 m s~^) to account for instrumental effects using the telluric lines. 

The final spectra were of high S/N reaching ^ 130 per pixel (880 A~^) in the H & K 
emission core, 400 per pixel (2700 A~^) in the pseudo-continuum near 3950 A, and about 
3 times higher near the Ca II IRT. Spectra with comparable S/N were taken of 61 Vir, a 
G5 V star known not to have close-in giant planets, plus the hot standard HR 5511 (SpT = 
AOV) for telluric line correction. Table 2 lists the program stars, including their magnitudes, 
exposure times, and typical S/N. 

All further processing and analysis were performed with standard IRAF (Image Reduc- 
tion and Analysis Facility) routineso Differential radial velocity corrections were applied to 
each stellar spectrum using IRAF's fxcor and rvcorrect routines. Representative spectra near 
the key stellar activity indicators are shown in Figured] for HD 189733. 



^Also see |http://www.c£ht. hawaii.edu/Instruments/Spectroscopy/Espadons/Espadons jesprit.html 

^IRAF is distributed by the National Optical Astronomy Observatories, which is operated by the As- 
sociation of Universities for Research in Astronomy, Inc. (AURA) under cooperative agreement with the 
National Science Foundation. 
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3. Measuring chromospheric activity 

The very strong Ca II H and K photospheric absorption lines suppress the local stellar 
continuum making it difficult to normalize each spectrum consistently. The normalization 
level was set at 0.3 of the flux at 3950 A centered on the H and K lines. Therefore the 
wavelengths were constant for all spectra of a given star, though they varied slightly from 
star to star due to variations in spectral type. The :^7-A spectral range was chosen to isolate 
the H and K reversals. This window is wide enough that a few photospheric absorption 
features appear to test for general stability. To normalize each sub-spectrum, the end points 
were set to 1 and fitted with a straight line. The mean Ca II K cores for four of the program 
stars are shown in Figure [2] with those for HD 179949, HD 189733 and v And in Figures [3] 

-El 

The spectra were grouped by date and a nightly mean was computed for each of the 
lines. Other than the active HD 189733 (see Section [3l2i) . all stars observed had non- varying 
K emission at the < 0.001 level on average on a given night. This is a result of S/N variations 
and intranight (short time-scale) chromospheric activity. 

We used nightly residuals from the average stellar spectrum to measure the chromo- 
spheric activity within the reversals. Each residual spectrum had a broad, low-order curva- 
ture removed. The residuals of the normalized spectra (smoothed by 17 pixels) were used 
to compute the mean absolute deviation (MAD = N~^'L\datai — mean\ for N spectra), a 
measure of overall variability within the span of the observing run. The Ca II K MAD 
spectrum and the nightly residuals used to generate it for HD 179949, HD 189733 and v 
And are displayed in Figures [31 — O 

The analysis presented in this section consists of only Ca II K emission measurements 
for several reasons: 1) The broad, deep photospheric absorption of the Ca II K line allows the 
chromospheric emission to be seen at higher contrast as compared to Ha and the Ca II IRT 
where chromospheric emission merely ffils in the absorption core. 2) Extensive studies by the 
Mount Wilson group and Wright et al. (2004) allow us to isolate the average chromospheric 
emission (K') by correcting for the photospheric contribution to our measurements ((K') = 
{K) — {Kphot), see Section 3.3.2 of Shkolnik et al. 2005a for more details.) This makes for a 
more accurate comparison of the stars in the sample whose spectral types vary from F7 to 
Kl. 3) Previous CFHT/Gecko spectra consisted of only a single order containing Ca II H & 
K and it is useful to make comparisons between the data sets. 4) Lastly, there are no telluric 
features or blended lines to contaminate the spectra, as is the case for the other indicators, 
which are discussed in Section HI 
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3.1. HD 179949 & v And: Evidence of the on/off nature of SPI 

When monitoring chromospheric emission, steUar activity may be modulated by the 
star's rotation, planetary motion in the case of SPI, or a combination of both. The orbital 
periods of the planets are well known and uniquely established by the PRV and transit 
discovery methods, but the rotation periods of the stars are much harder to determine in 
part due to stellar differential rotation. For studies of SPI, differentiating between rotational 
and orbital modulation of the chromospheric emission is key. 

In Shkolnik et al. (2005a) we presented evidence of planet-induced heating on HD 179949. 
The effect lasted for over a year and peaked only once per orbit, suggesting a magnetic 
interaction. In the simplest configuration, a magnetic interaction would occur near the sub- 
planetary point, when the planet is in front of the star relative to the line-of-sight, which 
defines orbital phase (porb = 0. Reproduced in Figure [6|, we fitted a truncated, best-fit spot 
model to our 2001 and 2002 data with P = P^rb = 3.092 d corresponding to the change in 
projected area of a bright spot on the stellar surface before being occulted by the stellar 
limb. The fit to the 2001 and 2002 data peaks at (porb = 0.83 ±0.04 with an amplitude of 
0.027. We over-plot new data from 2005 which is fit remarkably well by the same model 
with only an insignificantly small relative phase shift of -0.07. 

This phase lead may help identify the nature of the interaction. For example, the offset 
from the sub-planetary point of a starspot or group of starspots can be a characteristic 
effect of tidal friction, magnetic drag or reconnection with off-center stellar magnetic field 
lines. For further discussion on such mechanisms, see papers by Gu et al. (2005), Preusse et 
al. (2006) and Mclvor et al. (2006). In any case, the phasing, amplitude and period of the 
activity have persisted for over 4 years. 

Ca II data acquired in 2003 and 2006 of HD 179949 do not phase with the planet's orbit 
(Figure [7]), but both phase well with a 7-day period, likely the rotation period of the star. 
In Figure [S] we fit data from each year separately with a rotation curve because the effects 
of differential rotation and the appearance and disappearance of new spots over the three 
years would produce variations in phase, amplitude and period in the observed modulation. 
Note that the amplitude of the rotational activity in only 0.6 of that of induced by SPI. 
Indirect indications of the rotation rate of HD 179949 imply P^ot ~ 9 days and are presented 
in Shkolnik et al. (2003) and Saar et al. (2004). Wolf & Harmanec (2004) weakly detect 
(1.5 a) a photometric rotation period for HD 179949 of 7.07 d with an amplitude of only 
0.008 mag. While more photometry is needed to determine a rotation period conclusively, 
the modulated Ca II emission of this star in both 2003 and 2006 strongly suggests a rotation 
period of 7 days. 
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Similarly, previous Ca II data of v And indicated possible SPI (Figure 8 of Shkolnik 
et al. 2005a), yet our September 2005 data appears to vary with the rotation. Again, the 
rotation period is not well known. Henry et al. (2000) quotes both 11 and 19 days, with a 
probable 11.6-day period from the (Shk) index. We plot the 2005 data against the 4.6-d 
orbital in Figure [9] and an 11.6-day rotation period in Figure [TOl Unlike data from 2002 and 
2003, the 2005 data phase much better with Prot=^^-Q d than with a planetary orbit. 

This on/off characteristic of SPI observed in the HD 179949 and v And systems is 
predicted by the models of Cranmer Sz Saar (2007). They model the Ca II H & K light curve 
of a sun-like star with a hot Jupiter interacting with the field geometry at various stages 
of the empirically derived solar magnetic field at annual steps of the 11-year solar cycle. 
They conclude that due to the complex nature of the multipole fields, the Ca II K light 
curves due to SPI do not repeat exactly from orbit to orbit, and at times the planet-induced 
enhancement may disappear altogether leaving only rotationally modulated emission. This 
may explain the 2003 and 2006 disappearance of the strong orbital modulation seen in 2001, 
2002 and 2005 for HD 179949. Their models also show that for sparsely sampled data, the 
apparent phase shift between the peak Ca II emission and the sub-planetary point may fall 
between 0.2 and +0.2 (or ± 72°), consistent with the -0.17 phase shift we detected repeatedly 
for HD 179949. 



3.2. HD 189733: The active host of a massive planet 

We reported in Shkolnik et al. 2005a that chromospheric variability of the active, young 
star Ceti ((K') =0.815), for which the presence of a hot Jupiter is not ruled out, and 
HD 73256, known to host a hot Jupiter (MpSim=1.85 Mj, a=0.037 AU, Udry et al. 2003, 
(K')=0.899) was modulated by stellar rotation with additional variability or flaring poten- 
tially induced by a hot Jupiter. We find a similar effect on HD 189733, a generally more 
active star, for which the average Prot is well known (11.73±0.07 days, CroU et al. 2007). 
This star has relatively strong K emission with (K') =1.231 and large intranight variability 
as shown in Figure [11] which varies on time-scales at least as short as the length of the indi- 
vidual exposures (30 minutes). The average emission from night to night clearly varies with 
the star's rotation though with a lower amplitude (0.011 A) as compared to HD 73256 (0.045 
A), likely because the star is intrinsically more active, with a larger percentage of its surface 
covered in spots. This makes it difficult to extract any magnetic and/or tidal contribution 
to HD 189733's chromospheric emission by the planet as no correlation is seen between the 
planet's orbit and the residuals to the rotational modulation. However, Figure [12] shows 
how the integrated mean absolute deviation of the K line residuals from the global mean per 
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night (MADK) vary with orbital phase. Though we only have four nights of observations 
spanning 1.4 orbits, there is a clear increase in very-short-term (<30 minutes) activity at 
4>orb ~ 0.8. Remarkably, this is the same phase at which SPI peaks for HD 179949 and r Boo 
as measured both spectroscopically and photometrically. (See discussion below and Walker 
et al. 2006, 2007) 

3.3. r Boo: SPI on a tidally-locked star 

The star with the shortest rotation period in our sample is r Boo. It has the largest 
f sim (= 14.8 m s~^. Gray 1982) and is believed to be in synchronous rotation with its tightly- 
orbiting massive planet {Prot = 3.2 ± 0.5 d, Henry et al. 2000, Porb = 3.31250 d, MpSim=4.4 
Mj, Butler et al. 1997). We observed a small but significant night-to-night modulation in 
the H & K emission of r Boo during the first 3 years of observations with no obvious phasing 
with the planet's position. (Data from 2001, 2002 and 2003 are in Figure 7 of Shkolnik et 
al. 2005a.) Due to the tidal- locking of the star with the planet, we must depend on consistent 
orbital phasing of any modulated activity to disentangle SPI from rotational modulation for 
T Boo. 

Walker et al. (2007) presented light curves of r Boo taken in 2004 and 2005 observed 
in broadband optical light by the MOST space telescope. In the first year, they observed a 
significant photometric signal close to the planet's orbital frequency. In the second year, there 
was no signal of similar strength but a clear correlation with the MAD of the photometry 
with orbital phase. They showed that when phased to the planet's orbital period, the active 
region precedes the sub-planetary point by 68°, very close to the phase lead we observe in 
the enhanced activity on HD 179949 and in the MAD of HD 189733's Ca II K emission. 
Though synchronisity with the exact planet's position is not obvious in r Boo's 2001—2003 
Ca II data. Walker et al. (see their Figure 6) show that the MAD of these data during the 
photometrically active phase range, centered on (port =0.8, is twice as high than outside of 
it. Though of smaller amplitude (and larger error bars) than in previous years, our new 2006 
Ca II data (Figure [T3|l R may also show a weak enhancement between (f)orb =0.7 and 1.2. If 
this is indeed the case, it implies that an active region leading the sub-planetary point, has 
persisted on r Boo for at least 5 years, equivalent to ~550 planetary orbits. 



^Of the 9 June 2006 nights on which t Boo was observed, three were taken in 'spectropolarimetry' mode, 
those at orbital phases of 0.19, 0.49 and 0.82. It is interesting to note that on the two nights with high K 
emission, Catala et al. (2007) detect a clear Stokes V signature while the third night, (0orb=O.49) has both 
low K emission and no Stoke V signal. 
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3.4. Night-to-night activity correlates with planet's magnetic moment 

Sanchez-Lavega (2004) looked at the internal structure and the convective motions of 
giant extrasolar planets in order to calculate their dynamo-generated surface magnetism. 
Given the same angular frequency (which is a reasonable approximation for the short-period 
planets in question), the magnetic dipole moment, and hence the magnetospheric strength, 
increases with planetary mass. This is observed for the magnetized planets in our own so- 
lar system, where the magnetic moment grows proportionally with the mass of the planet 
(Stevens 2005), and more specifically, with the planet's angular momentum (L oc MpRpP~J:^^, 
Arge et al. 1995). Since only lower limits exist for the masses of most hot Jupiters and at 
such small semi- major axes they should be tidally locked {Pp^rot = Port), we plot MpSini/Porb 
against (MADK), the average of the integrated MAD of the K line residuals per observing 
run, in Figure [TH Though we are able to include only one additional active point (HD 
189733) to the original plot of Shkolnik et al. 2005a, we continue to see an intriguing corre- 
lation between the planet's magnetic moment and the night-to-night chromospheric activity 
on its star. Of our sample, r Boo has the most massive planet and yet falls well below the 
correlation. This is consistent with the proposed Alfven wave model where the near zero 
relative motion due to the tidal locking of both the star and the planet (P*,roi = Pp,rot = 
Porb) produces minimal SPI because of the weak Alfven waves generated as the planet passes 
through the stellar magnetosphere, thereby transporting little excess energy to the stellar 
surface along the magnetic field lines (Gu et al. 2005). If this correlation between short-term 
activity and planetary magnetic moment holds for more hot Jupiter systems engaging in 
SPI, this could provide an empirical tool with which to estimate the strength of extrasolar 
planetary magnetic fields. 



3.5. Long-term stellar activity cycles 

With CFHT data spanning 5 years, we can compare the long-term variations in the 
chromospheric level of HD 179949, r Boo, v And and HD 209458. We measure Ca II K 
emission strength (K) by integrating across the normalized K cores bounded by the Kl 
features (Montes et al. 1994) and plot their average for each observing run in Figure [13 The 
error bars represent the MAD value for each observing season. There does not appear to 
be any correlation between the mean chromospheric activity and the level of night-to-night 
modulation, be it due to SPI or stellar rotation. Though 5-year baseline is a good start to 
tracking the intrinsic stellar activity cycles of these stars, limiting their periods to >10 years, 
the variability from run to run may also be due to active regions on the visible disk of the 
star. We require more frequent monitoring over several more years to firmly say anything 
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more about the activity cycle of any individual program star. 

4. Correlating with other activity indicators 

A similar analysis to the Ca II K line described above was performed for the Ca II H, IRT 
line at 8662A and Ha though with the normalization points set at 0.7 of the local continuum 
for the latter two. This level is within the photospheric damping wings of the line, focusing 
the analysis on the chromospheric core and excluding some blended and telluric lines. For 
both HD 179949 and HD 189733, there is a strong correlation between the residuals of Ca 
II K with those of the Ca II H and 8662A lines (Figure [TBI) as expected given the common 
upper energy level of their transitions. However, a poorer correlation exists with the Ha lines. 
Though Ha is often demonstrated to be just as good a tracer of chromospheric activity as Ca 
II, a recent analysis by Cincunegui et al. (2007) has shown that when comparing a sample 
of stars, the correlation between Ha and Ca II is the result of the correlation of each line 
with spectral type rather than with stellar activity. When comparing the variability of the 
lines for an individual star, there is no consistent correlation. This is likely the effect of 
the differing underlying formation physics between them. (Soderblom et al. 1993, Cram & 
Mullan 1985). We leave the relative energy emitted in the lines and their implications for 
SPI models for a later paper. 

The He I D3 line (a blended triplet) at 5876 A correlates well with plage regions on the 
solar surface (Landman 1981) as well as with the Ca II H&K emission (Saar et al. 1997). 
It forms in the upper chromosphere and is thought to be back-heated by the stellar corona 
giving us a unique optical view of this hot plasma. The absence of He I absorption in non- 
magnetic regions on the sun and in other inactive stars indicates that He I has no basal 
(acoustically heated) flux level, unlike the other activity indicators in the visible spectrum, 
and is therefore purely a signature of magnetic activity. 

We show the spectral region near the He I lines of our program stars in Figure [T71 
The line is blended with the weak lines of Fe I 5876.30, Cr I 5876.55, and unidentified lines 
at 5875.76 and 5875.14 A, both blended with telluric H2O (Moore et al. 1966). The im- 
perfect removal of the telluric lines from our spectra left residuals at the level of <0.003 
of the nearby continuum. This made it difficult to analyze night-to-night variations in the 
stars that exhibit relatively strong variability in the Ca II lines, though the mere presence 
of the line indicates magnetic heating source, with strong absorption implying great activ- 
ity. HD 179949, HD 189733 and r Boo all have clear He I absorption while HD 209458, 
HD 149143, HD 217107 as well as our standard star 61 Vir have only weak, if any, absorp- 
tion. Though night-to-night variability is difficult to quantify in HD 179949 and HD 189733, 
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there is a clear increase in He I D3 absorption on the night that each of the two stars displays 
its maximum Ca II emission, which do not occur on the same night |f| 

We measured the average He I EW by deblending it with the contaminated lines, a 
technique particularly difficult for the rapid rotator r Boo. The values are listed in Table [1] 
It is interesting to note that the He I EW does not follow the power-law correlation with Ca II 
H & K emission observed by Saar et al. (1997) for G and K dwarfs, but does display a strong 
correlation with the short-term activity metric (MADK) plotted in Figure [TSl where stars 
with more night-to-night activity have stronger He I absorption. The relationship between 
the He I absorption and stellar rotation period observed by Saar et al. is also not obvious 
in our small sample. It therefore remains possible that the strength of the He I D3 line may 
predict whether a planet-bearing star will have night-to-night variability. When more stars 
with hot Jupiters are discovered this may come in handy in helping to decide which systems 
should be studied further with intensive time-series observations. 



5. Summary 

We have observed 7 stars with hot Jupiters using CFHT's echelle spectrograph ES- 
PaDOnS to search for night-to-night modulation of the Ca II emission for evidence of SPI. 
Four of these have been observed in our previous studies of time-varying Ca II H & K: HD 
179949, V And, HD 209458, and r Boo, and we have added three new targets: HD 189733, 
HD 217107 and HD 149143. 

For our prime target, HD 179949, we now have a total of six observing runs spanning 5 
years. During four runs (Aug 2001, July 2002, Aug 2002 and Sept 2005) the Ca II emission 
varied with the orbital period of 3.092 days, with consistent amplitude and peak phase 
indicative of a magnetic interaction between the star and planet. The peak activity on HD 
179949 in these epochs occurs at (port ~ 0.8, leading the sub-planetary longitude by some 
70°. Interestingly, this same phase shift is observed in the MAD of the Ca II K residuals of 
both r Boo and HD 189733. The phase lead can provide information on the field geometries 
(i.e. Parker spiral) and the nature of the effect such as tidal friction, magnetic drag or 
reconnection with off-center magnetic fields. 

HD 179949 data from the other two runs (Sept 2003 and June 2006) clearly vary with 



■^We cannot measure the He I EW of HD 179949 from the September 2005 data since poor weather 
prevented us from observing a telluric standard. Similarly, a telluric standard was not observed for the 
spectropolarimetric study of v And. 
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the rotation period of 7 days. A similar effect is seen on v And where one of four epochs 
appears to be modulated by rotation rather than the planet's motion. This on/off hehayior 
has been modeled by Cranmer & Saar (2007) to be an effect of magnetic reconnection with 
the stellar field as it varies with the star's long-term activity cycle. 

We present the expected correlations between the variability observed in the Ca II K 
line with Ca II H and IRT 8662, and a weaker correlation with Ha. Though we could not 
accurately measure the variability in the upper- chromosphere line He 1 D3, we show that it 
has the potential to flag stars which might be active in Ca II K in on a night-to-night time 
scale. 

To date wc have observed 13 stars with hot Jupitcrs at CFHT and VLT, of which 
5 appear to be actively engaging in SPl: HD 179949, v And, HD 189733, HD 73256 and 
r Boo. The activity as measured by the mean absolute deviation over a run on the first 
four of these stars correlates well with Mpsini / Pgrb, a value proportional to the planet's 
magnetic moment, and thus with the hot Jupiter's magnetic field strength. Because of their 
small separation (< 0.075 AU), a hot Jupiter lies within the Alfven radius of its host star, 
allowing a direct magnetic interaction with the stellar surface. Although this correlation is 
tentative, short-term chromospheric variability may be our first probe of extrasolar planetary 
magnetospheres . 

We are very grateful to Claude Catala and John Barnes for contributing several r Boo 
and V And spectra, to Jean-Francois Donati for making Libre Esprit available to CFHT 
users, to the CFHT staff for their excellent support of this program, and to Benjamin Brown 
for useful discussions on planetary dynamos. Research funding from the NASA Postdoctoral 
Program (formerly the NRC Research Associateship) for E.S. and the National Research 
Council of Canada (D.A.B.) are gratefully acknowledged. 
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Fig. 1. — Selected regions of a normalized spectrum of HD 189733 identifying key stellar 
activity diagnostics. 
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Fig. 2. — The mean normalized Ca II K cores for four of the six program stars. That of HD 
179949, HD 189733 and v And are shown in Figures [3l H] and [5l respectively. 




Fig. 3. — Top: The mean normalized Ca II K emission of HD 179949 on the three nights 
observed in September 2005 and the four nights in June 2006. Middle: The residuals relative 
to their respective means. Bottom: The mean absolute deviation (MAD) of the residuals. 
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Fig. 4. — As in Figure [3] for HD 189733 observed in June 2006. 
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Fig. 5. — As in Figure [3] for v And observed in September 2005. 
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Fig. 6. — Integrated flux of tlie K-line residuals from a normalized mean spectrum of 
HD 179949 as a function of orbital phase for the 2001 and 2002 data (open circles) pub- 
hshed in Shkolnik et al. (2003) and 2005 data (filled circles). The grey hne is a best- fit spot 
model to the earher data whose thickness reflects the error in the phase shift. The black hne 
is the same fit slightly shifted in phase by -0.07 to better fit the 2005 data. This small shift 
relative to the earlier data is not signficant. Error bars in the integrated residual K flux are 
± the intranight residual RMS. Note that the data points are repeated for two cycles. 
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Fig. 7. — Integrated flux of tlie K-line residuals from a normalized mean spectrum of 
HD 179949 for 2003 and 2006 data plotted on the 3.092-day orbital period with the SPI 
spot model from Figure [H] over-plotted. Error bars are ± the intranight residual RMS. 
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Fig. 8. — Integrated flux of tfie K-line residuals from a normalized mean spectrum of 
HD 179949 for 2003 and 2006 data plotted on a 7-day rotation period with phases rela- 
tive to the first night of each run. The points are vertically shifted such that the minimum 
of each curves is zero. Error bars are ± the intranight residual RMS. The curves are best- fit 
spot models to the two data sets. 
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Fig. 9. — Integrated flux of the K-line residuals from a normalized mean spectrum of v And 
for September 2005 data plotted on the 4.6-day orbital period. Error bars are within the 
size of the of the points. The curve is the best-fit spot model to our 2002 and 2003 data. 
(See Figure 8 of Shkolnik at al. 2005.) 




Fig. 10. — Integrated flux of the K-line residuals from a normalized mean spectrum of v And 
for September 2005 data plotted on a 12-day rotation period. Error bars are within the size 
of the of the points. The curve is a best-flt spot model. 
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Fig. 11. — Integrated flux of tlie K-line residuals from a normalized mean spectrum of 
HD 189733 as a function of a 11.7-day rotational phase (Moutou ct al. 2007, Croll et al. 2007). 
The open circles arc residuals of individual exposures relative to a global mean, while the 
solid circles use only the mean of each night. Error bars represent twice the measurement 
error of a given exposure as measured by the residuals outside of the Ca II K core. The solid 
curve is a best-fit sinusoid tracing the rotation of the star. 
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Fig. 12. — Mean absolute deviation (MAD) of the nightly K-line residuals shown in Figure [TT] 
of HD 189733 as a function of the 2.2-day orbital period. The error bars represent the 
integrated MAD immediately outside the Ca II emission core and reflect the S/N obtained 
for each night. Errors in phase are less than the size of the points. 
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Fig. 13. — Integrated flux of the K-line residuals from a normalized mean spectrum of r Boo 
as a function of orbital phase with a best-flt spot model. Note that since Prot — Porh foi' this 
system, 0,.ot ~ (borb- The three observations taken in 'spectropolarimetry' mode are at (prot 
— 0.19, 0.49 and 0.82. Error bars are ± the intranight residual RMS. 
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Fig. 14. — The ratio of the minimum planetary mass (in Jupiter masses) to the orbital period 
(in days) plotted against the average MAD of the K-line per observing run for all 13 stars 
observed. The x-axis quantifies the planet's magnetic moment assuming tidal locking, such 
that Prot = Porb- The fiUed-in circles are of stars which exhibit significant night-to-night 
variability in the Ca II K fine: HD 73526, v And, HD 179949 and HD 189733 (this work). 
The tidally-locked system of r Boo does not follow the correlation traced out by the others. 
The error bars are one-sided due to the positive contribution of integrated MAD immediately 
outside the Ca II emission core and reflect the S/N obtained for each target. 
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Fig. 15. — The average integrated Ca II K emission of HD 179949, r Boo, v And and HD 
209458 for each of the five or six observing runs. The error bars are the MAD values for 
each observing season. 
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Fig. 16. — Integrated flux in the Ca II K residuals plotted against the residual flux in Ca II 
H, the 8662A hne and Ha for HD 179949 and HD 189733. The number in the corner is the 
correlation coefficient of the best-fit line. 
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Fig. 17. — Normalized spectra in the region of the He I line (5876 A) for 5 program stars 
plus the standard, 61 Vir. The vertical line corresponds to the location of the line. The 
same spectral region for HD 189733 is shown in Figure [H The He I cannot be measured for 
V And since no telluric standard was observed on those nights. 
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Fig. 18. — Average MAD of the K emission (MADK) for the 2006 program stars as a function 
of the He I equivalent width. The error bars reflect the difficulty in measuring the EW due 
to line blends, especially for the rapidly rotating r Boo. 



